Images of Bursting Sources of High-Energy Cosmic Rays. I: Effects of 

Magnetic Fields 
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ABSTRACT 

It has recently been shown that the highest energy cosmic rays (CRs) may originate 
in the same cosmological objects producing 7-ray bursts. This model requires the 
presence of intergalactic magnetic fields (IGMF) to delay the arrival times of ~ 10 20 
eV CRs by 50 years or longer relative to the 7-rays, of an amplitude that is consistent 
with other observational constraints. Sources of CRs coming from individual bursts 
should be resolved with the planned "Auger" experiment, with as many as hundreds 
of CRs for the brightest sources. We analyze here the apparent angular and energy 
distribution of CRs from bright sources below the pion production threshold (in the 
energy range 10 19 eV < E < 4 x 10 19 eV) expected in this model. This observable 
distribution depends on the structure of the IGMF: the apparent spectral width AE 
is small, AE/E <^ 1%, if the intergalactic field correlation length A is much larger 
than IMpc, and large, AE/E = 0.3, in the opposite limit A <C IMpc. The apparent 
angular size is also larger for smaller A. If the sources of CRs we predict are found, 
they will corroborate the bursting model and they will provide us with a technique to 
investigate the structure of the IGMF. 



Subject headings: cosmic rays — gamma rays: bursts — magnetic fields 



1. Introduction 



The origin of cosmic rays (CRs) with energy E > 10 19 eV is unknown. Most of the sources 
of cosmic rays that have been proposed have difficulties in accelerating CRs up to the highest 
observed energies (e.g., Cronin 1992 , Elbert &; Sommers 1995 , Sigl, Schramm fc Bhattacharjeaj 



1994 ). Recently, a new model has been proposed where the same astrophysical objects responsible 



for 7-ray bursts (GRBs) also produce the highest energy CRs ( Waxman 1995a , Vietri 1995] , 
Milgrom & Usov 1995| ). There is increasing evidence that GRBs and the highest energy CRs are 



of cosmological origin QFishman fc Meegan 199*5 , Bird et al. 1994, Yoshida et al. 1995, Waxman 



1995b). In this case, the observational characteristics of GRBs impose strong constraints on the 



physical conditions in the 7-ray emitting region, which make it an ideal site for accelerating 
protons up to the highest observed CR energies. In addition, the average rate (over volume and 
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time) at which energy is emitted as 7-rays by GRBs and in CRs above 10 19 eV in the cosmological 
scenario is remarkably comparable ( Waxman 1995aj ,b). 



The energy loss of CRs with E ^ 10 20 eV, due to interaction with the microwave background, 
imply that they must originate at distances < 100 Mpc (e.g., Pronin 1992| ). To reconcile the 
detection of two CRs with E £ 10 20 eV over a ~ 5yr period ( [Bird et al. T994] , |Yoshida et al. 1995| ) 



with the expected rate of nearby GRBs (~ 1 per 50yr in the field of view of the CR experiments 



out to lOOMpc; e.g., |Cohen fc Piran 1995|) , a dispersion in arrival times > 50yr for CRs produced 
in a single burst needs to be invoked. Such dispersion may result from deflections of CR protons 
by the intergalactic magnetic field (IGMF) ( Waxman 1995a , Waxman fe Coppi 1996j ). The 



deflection angle for a proton propagating a distance D in a magnetic field B with correlation 
length A is 9 S ~ 0.°025, (D/X) 1 / 2 (A/lOMpc) (B/l()- u G) (-E/10 20 eV)-\ and the induced time 
delay is t(E) ~ DO 2 /Ac ~ 200yr(Z?/100Mpc) 2 (A/lOMpc) {B/lt)- 11 G) 2 (£/10 20 eV)- 2 (where the 
exact dependence of r on 9 S depends on the scale of the magnetic field; see below). The required 
magnetic field is consistent with observational limits ( [Kronberg 1994 , Vallee 1990[ ). 



In a preceding paper ( Miralda-Escudc fc Waxman 199!f ), we calculated the number of CR 



sources with different flux at different energies in the model of "Cosmic Ray Bursts" . The energy 
dependent delay of the CR arrival times results in individual sources having relatively narrow 
observed spectra, since at any given time only those CRs having a fixed time delay are observed. 
It was shown that in the GRB model most of the CRs with E ;> 2 x 10 20 eV should come only from 
a few sources, since the short time delays implied by the small distance to the sources allow only 
a few bursts to be seen at any given time. At lower energies many more sources are observable 
due to both the larger distances over which the CRs can propagate and the much longer time 
delays. Therefore, while at the highest energies bursting sources should be identified from only 
a small number of CRs from their coincident directions, many more CRs need to be detected at 
low energies to identify sources. Recently, the AGASA experiment reported the presence of 3 
pairs of CRs with angular separations (within each pair) < 2.5°, consistent with the measurement 
error, among a total of 36 CRs with £>4x 10 19 eV ( |Hayashida et al. 1996| ). The two highest 
energy AGASA events were in these pairs. Given the total solid angle observed by the experiment, 
~ 2-7rsr, the probability to have found 3 pairs by chance is ~ 3%; and, given that three pairs 
were found, the probability that the two highest energy CRs are among the three pairs by chance 
is 2.4%. Therefore, this observation favors the bursting source model, although more data are 
needed to confirm it. 

The detailed spectral and angular distribution of CRs from an individual source, expected in 
a bursting source model, depend on the structure of the IGMF and on the details of the energy 
loss processes. This shall be examined in this Letter, confining the discussion to sources observed 
at energies below the pion production threshold, £^4x 10 19 eV. Sources at higher energy will be 
treated in a subsequent publication. 
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2. The Propagation of CRs through the Magnetic Field 

2.1. Large Scale Fields 

Let us first consider the case where the typical displacement of the CRs is much smaller 
than the correlation length, A 3> D9 S (D, E). In this case, all the CRs that arrive at the observer 
are essentially deflected by the same magnetic field structures. Therefore, if the CR energy loss 
during propagation was perfectly continuous and deterministic, all CRs with a fixed observed 
energy would reach the observer with exactly the same direction and time delay. At a fixed 
time, the source would appear monoenergetic and point-like, with the CR energy decreasing with 
time as i _1//2 , and the deflection from the true source position increasing as i 1//2 . However, the 
fact that CRs suffer random energy losses as they propagate, owing to the production of pions 
and pairs in interactions with the microwave background, implies that in reality the image 
should be blurred both in arrival direction and energy. Due to the random energy loss, at any 
distance from the observer there is some finite spread in the energies of CRs that are observed 
with a given fixed energy. Since the deflection due to the magnetic field is everywhere inversely 
proportional to the CR energy, the fractional variation of the total deflection, 59(E, D)/9 S (E, D), 
and of the time delay, 5t(E, D)/r(E,D), are similar to the characteristic fractional energy spread 
along the path, ue{E,D)/E. The spread in energy of CRs observed at a given time is also 
AE/E ~ St/r ~ a E {E, D)/E. 

The brightest sources of CRs with E ^ Ax 10 19 eV should be very nearby. The typical distance 
D m {E) to the brightest source observed over an energy range AE around E, with AE/E ~ 1, is 
the radius of a sphere within which the average time between bursts is equal to the characteristic 
time delay t[E, D m (E)]; i.e., D m is determined by A^D^yriE, D m )/3 = 1 where v is the burst 
rate per unit volume. Thus, 

DmiE) " 30 ( lO-«Mpc- 3 )" 1/5 ^ 9/5MpC ' (1) 

where E = 10 19 -EigeV and To = t(E = 10 20 eV, D = lOOMpc). The sources of the events that have 
been detected above 10 20 eV must lie within lOOMpc, requiring vtq > 10~ 6 Mpc~ 3 . From eq. (||) we 
see that at lower energies the brightest sources will be even closer. These CRs from bright sources 
will suffer energy loss only by pair production, because at E < 5 x 10 19 eV the mean-free-path for 
pion production interaction (in which the fractional energy loss is ~ 10%) is larger than lGpc. 
Furthermore, the energy loss due to pair production over lOOMpc propagation is only ~ 5%. 

The fractional energy loss in a pair production interaction is ~ m e /m p . The average 
total energy loss after TV" interactions is 5E/E ~ Nm e /m p , and the dispersion in the energy 
loss is oe/E ~ y/~N m e /m p = (m e /m p 5E / E) 1 / 2 . The spectral width and apparent angular 
size of bright sources at E < 4 x 10 19 eV, for which SE/E g 5%, are therefore very small, 
AE/E ~ 59 /9 S ~ cte/E 1%. On the other hand, CRs of much higher energy can lose a 
substantial fraction of their energy in a single pion production, resulting in a large dispersion in 
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the energy loss for a propagation distance of a few tens of megaparsecs. This implies that the 
width of the energy spectrum and the angular size of the image should be much wider at high 
energies. Detailed predictions of the spectral shape due to the random energy loss caused by pion 
production will be presented in a subsequent publication. 



Let us now consider the case where the typical displacement of the CRs is much larger than 
the correlation length, A <C D9 S (D,E). In this case, the deflection of different CRs arriving at 
the observer are essentially independent. Even in the absence of any energy loss there are many 
paths from the source to the observer for CRs of fixed energy E that are emitted from the source 
at an angle 6 <J 9 S relative to the source-observer line of sight. Along each of the paths, CRs are 
deflected by independent magnetic field structures. Thus, the source angular size would be of 
order 6 S and the spread in arrival times would be comparable to the characteristic delay r, leading 
to AE/E ~ 1 even when there are no random energy losses. 

To allow for an analytical treatment of the case 6 S D 3> A, we consider CRs below the pion 
production threshold and neglect their energy losses due to e ± pair production (which is, as shown 
above, a good approximation for bright sources of CRs with E <^ 4 x 10 19 eV). Since each CR 
arriving at the observer follows a different path, different CRs sample different "realizations" of 
the IGMF. The distribution of observed deflections and time delays for an observer at a distance 
r from the source is therefore given by the probability density P(<p, t; r, E) for a CR propagating 
to a distance r to encounter a magnetic field structure that would cause a time delay t = T — cr 
(where T is the time since the CR was emitted) and deflection angle tp relative to the direction to 
the source. 

Let us divide r in small intervals Ar such that A < Ar <C r. The deflections in different 
intervals are independent, since Ar is larger than the field correlation length. The evolution of 
the probability density P may then be described as a Markov process that changes P at each 
interval Ar. Since the number of independent deflections, r/ Ar, is very large, the relative change 
in P over Ar is small, and we may approximately describe its evolution by a differential equation 
obtained by formally taking the limit Ar — ► 0. 

In what follows we replace the variable if with £ = tp 2 , P(C, t; r, E) = P((p, t; r, E)dip/dQ. Since 
d( = cpdip/2, P(£,t;r, E) is the probability per unit solid angle. In the absence of magnetic field, 
as a CR moves to a distance r + Ar its angle <p decreases by Aip = —ipAr/r, i.e., A( = — 2£ Ar/r, 
and its time delay increases by At = (Ar/2c. Thus, in the absence of magnetic deflection the 
evolution of P is determined by 



where the right-hand-side is due to the change in the volume element dC, dt for CRs moving from r 



2.2. 



Small Scale Fields 




(2) 
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to r + Ar. In the presence of an IGMF, a CR has a probability f(a, Ar)da to be scattered by an 
angle a to a + da while propagating a distance Ar. CRs moving at an angle ip after a scattering 
by an angle a may have been scattered from a direction ip' obeying ip' 2 = ip 2 + a 2 — 2(£>acosp, 
where p is the angle over a circle of radius a around the direction ip. The average over this circle 
of C — £ is a 2 , and the average of (£' — C) 2 is 2£a 2 (where we use the approximation a <C ip, valid 
for our case of many scatterings along the CR trajectory). Therefore, the difference between the 
number of CRs scattered into (p and the number scattered out of this direction over the interval 
Ar is < a 2 (Ar) > (dP/dC, + Cd 2 P/d( 2 ), where < a 2 (Ar) > denotes the average of a 2 over the 
distribution /(a, Ar). Combining this expression with (g), our final equation is 

dP 2P /2C A dP d 2 P (dP 

where ip =< a 2 (Ar) > /Ar. The effect of the magnetic field is therefore to introduce a 
"mean scattering rate," ip. If the IGMF is characterized by a single correlation length A, then 
ip =<a 2 (X) > /A = 2(eB/E) 2 X/3, where the 2/3 factor accounts for random field orientations. If 
the field fluctuates on a range of scales A, with some power spectrum dB 2 (X) = &\d\, then 

^ ~ U / dl ° gX * xX * = Ifr^ • (4) 

Equation (|3|) is analogous to the equation describing the propagation of photons in a scattering 
medium in the limit of infinite optical depth, and was derived for the latter case by Alcock & 
Hatchett (1978), who presented an analytical solution (their eqs. [21]-[24]). Of particular interest 
to us is the integration of P over £, giving the distribution of time delays. We can express this 
in terms of the energy distribution of CRs with fixed time delay, which is the most relevant 
observable quantity given the long time delays expected in the bursting model. The differential 
flux received from a burst, which produced N{E)dE CRs in the energy range E to E + dE, is: 

E 2 



F(E; D, t) = _ E 2 N(E) V (-l) n+1 n 2 exp 



-2nV 



E$(t,D) 



(5) 



where Eo(t,D) = De(2B 2 X/3ct) 1 ^ 2 . Current observations require the generation spectrum N(E) 
to satisfy N(E) oc E~ n with n ~ 2 ( |Waxman 1995b| ). Therefore, E 2 N(E) on the right hand 



side of (||) is approximately constant, and the differential spectrum of the source is given by 
F(E)/E 2 N(E), presented in Fi gure 1. For this spectrum, the ratio of the rms CR energy spread 
to the average energy is 30%. 

The form of the angular image of CRs is also predicted. The angular distribution of CRs 
of a fixed energy averaged over all time delays is a gaussian, P((;E) oc exp[— (/9 2 (E)] with 
9 S {E) = e{2DB 2 \/9) 1 / 2 /E. The angular distribution of CRs at a fixed time delay averaged over 
all energies, which is more interesting given the long time delays expected, is approximately given 
by P((;t) oc exp{-0.58[C/6» 2 (^)] 2 }, where E is the average energy of CRs observed at a fixed time 
delay t. The joint energy and angular distribution of CRs at a fixed time delay, which may be 
observable, can also be predicted in detail from the solution of Alcock & Hatchett (1978). 
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3. Discussion 

We have shown that in a bursting source model the spectral shape and angular size of bright 
CR sources below the pion production threshold, E 4 x 10 19 eV, are sensitive to the ratio of 
the typical displacement of CRs, 9 S D, and the IGMF correlation length A. The spectral width 
AE and angular size A9 are small, AE/E ~ A9/9 S <^ 1%, in the limit 6 S D/X <C 1, and large, 
AE/E ~ A9/9 S ~ 1, in the opposite limit 9 s D/\ 3> 1. The spectral shape of the sources in the 
limit 9 S D/X ^> 1 is given by eq. (||) and shown in Figure 1. In this limit, the number of CRs 
observed at a fixed time delay with deflection <p to </? + dip relative to the direction to the source 
is approximately proportional to 99exp{— 0.58[(p 2 /9 2 (E)] 2 }dip, where E is the average energy of 
CRs observed at the specified time delay. If magnetic fields were strong in intermediate scales, 
9 S D/X ~ 1, a single source could have several images with angular separation ~ 9 S , which should 
be magnified and distorted. As in the case of gravitational lensing, an odd number of images 
should be expected with highly magnified images appearing in merging pairs near caustics of the 
magnetic deflection mapping. Each image could have a narrow energy spectrum AE/E <C 1, but 
usually the magnetic field on much smaller scales would be important in widening each image 
both in angle and energy. The average energy of CRs in different images would differ by AE ~ E. 

The condition 9 S D = 1 may be written as 

B 10-10 (-JLY V2 El9 GMpc" 1 / 2 . (6) 



V30Mpc 

Thus, measuring the spectral width of bright CR sources would allow to constrain the IGMF 
parameters provided the distance scale to the bright sources is estimated. In a bursting source 
model, the typical distance to the sources is proportional to (vtq)~ 1 /^ , where v is the bursting 
rate per unit volume and tq = t(E = 10 20 eV, D = lOOMpc) the characteristic time delay (cf. eq. 
]!]]). Current observations require vtq > 10~ 6 Mpc, which sets an upper limit to the distance of 

2/5 

the brightest sources, D m {E) <^ 30-E^g Mpc. Future CR experiments would allow to obtain better 
constraints on vtq by analyzing the number of sources at different energies ( Miralda-Escude~fe| 



Waxman 1996| ). Furthermore, if bursts originate in galaxies (as it seems likely if they are GRBs), 



their parent galaxies should be identified given that they are nearby, and the distance to each 
burst would then be known. For the GRB model, where v ~ 3 x 10 _8 Mpc _3 yr- 1 , the upper 
limit on BX 1 ! 2 from Faraday rotation measures, BX 1 / 2 < 1(T 9 G Mpc 1/2 (|Kronberg 1994], |Vallea; 



1990[) , sets an upper limit tq < 10 5 yr to the time delay and a lower limit to the distance scale 



2/5 

D m {E) p> 10-Eqg Mpc. Figure 2 shows the line 9 S D = A in the B — X plane together with the 
Faraday rotation upper limit BX 1/2 < 1CT 9 G Mpc 1/2 and the lower limit BX 1 / 2 > lCT n G Mpc 1/2 
required in the GRB model to allow tq > 50yr. The 9 S D = X line divides the allowed region in the 
B — X plane at A ~ IMpc. Thus, measuring the spectral width of bright CR sources would allow 
to determine if the field correlation length is much larger, much smaller, or comparable to IMpc. 



The proposed Auger experiment ( |Cronin 1992| , [Watson 1993| ) should be able to detect many 



CRs from individual sources at energies of a few times 10 eV. In fact, the number of CRs produced 
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by the brightest source at energy E and detected by a detector with area A and exposure time T is 
approximately N m (E) ~ fATN(> E)/4irD m (E) 2 T(D m , E), where iV(> £) is the number of CRs 
produced by the source above -E and / the fraction of time during which the source is within the field 
of view. Using eq. (0) and vN(> E) ~ 5 x lO^f^Mpc" yr 1 , as required by current observations 
([Waxman 1995b|) , we find iV m (£) ~ 10 3 J B 1 " 9 3/5 (^ro/10" 6 Mpc _3 )" 1 / 5 (/AT/10 4 km 2 yr). The Auger 
experiment will have 5000km 2 detectors and energy resolution of ~ 10%, giving an expected number 
of events detected from individual bright sources 10 2 — 10 3 (assuming a burst rate comparable 
to that of cosmological GRBs), which should easily allow to determine the spectral width of the 
sources. Furthermore, with an expected angular resolution of ~ 1°, the apparent angular size 
of the sources may be measured in the case 9 S D ^> A: the angular size of the brightest source 
at energy E is expected to be ~ 6 s [D m (E),E] ~ 0. 7S 19 4/5 (^r /10- 6 Mpc- 3 )- 1 /iO( ro /50y r )i/2. 

Measuring the angular size of the sources would allow to break the utq "degeneracy," and to 
estimate separately the characteristic time delay To oc -BA 1 / 2 and the burst rate v. 

Finally, we note that the deflections due to the magnetic field in the Milky Way will be 
superposed with those of the IGMF. While the additional time-delay caused by our galaxy is 
negligible, the deflection is not. For example, for a 3fiG field over a scale of 1 Kpc and a CR 
energy 2 x 10 19 eV, the deflection is 8°. This means that if the IGMF is present only on small 
scales, large coherent deflections could still be caused by the galactic field, which could produce 
elongated images with the CR energy being correlated with the observed deflection from the true 
source position, as well as multiple images of individual bursts. 

This research was partially supported by a W. M. Keck Foundation grant and NSF grant 
PHY95-13835 to the IAS. 
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Fig. 2. — The line 6 S D = X for a source at a distance D = 30Mpc observed at energy E ~ 10 eV 
(dot-dash line), shown with the Faraday rotation upper limit BX 1 / 2 < 10~ 9 G Mpc 1 / 2 (solid line), 
and with the lower limit BX 1 / 2 > 10 _11 G Mpc 1 / 2 required in the GRB model to allow tq > 50yr 
(dashed line). 



